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Introduction 

Aeroelasticity that involves strong coupling of fluids and structures is an important 
element in the design and optimization of aerospace vehicles. Large general purpose 
codes such as ASTROS[l] and NASTRAN[2] can compute the aeroelasticity of complex 
geometries using finite element (FEM) structures and linear aerodynamic equations. 
These codes use the AIC (aerodynamic influence coefficient) approach to couple fluids 
and structures. The concept of using AIC for aeroelastic computations was introduced 
several years ago [3]. The piston theory proposed for supersonic flows[4] introduced 
the AIC concept that became suitable with the finite element structural analysis 
approach. The early effort in developing the AIC approach compatible with FEM was 
presented by Appa et. al.[5 ] 


Several real world examples have shown that linear aerodynamic theory is not adequate 
for accurately predicting the aeroelasticity’ associated with complex flows. Modem 
supersonic transports that have highly swept wings can experience vortex- induced 
aeroelastic oscillations [6]. In order to predict this phenomenon, a capability to directly 
couple the high fidelity Euler/Navier-Stokes (ENS) equations with modal structures was 
needed[7]. As another example, the aeroelasticity of space planes that involve strong 
coupling of fluids, structures and controls is an important element in the design process 
An instability can occur soon after the space plane is separated from its carrier. The 
phenomenon is dominated by complex flows coupled with structural motions. From the 
results presented in Ref. 8 it is observed that the linear aerodynamic method was not 
adequate to completely understand the instability phenomenon which involved non-linear 
flows coupled with structural motions. High fidelity equations such as Euler/Navier- 
Stokes (ENS) for fluids directly coupled with finite elements (FE) for structures are 
needed for accurate aeroelastic computations in which these complex fluid/structure 
interactions exist. Use of high- fidelity equations involves additional complexities that 
result from numerics such as higher-order terms. Therefore, the coupling process is more 
elaborate when using high fidelity methods than it is for calculations using linear 
methods. 
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In recent years efforts have begun to improve the fidelity of AIC using non-linear CFD 
computations in place of linear aerodynamics. In Ref. 9 an approach to correct the AIC 
approach based on linear aerodynamic theory by using the steady state Navier-Stokes 
solutions is presented for cases with small deflections. Reference 10 presents use of the 
Euler equations to compute AIC using transpirent boundary conditions for cases with 
small deflections. In both approaches, AICs are computed without deforming CFD grids. 
For accurate computations of AIC which are valid for real world cases it is necessary to 
deform CFD grids. Efficient CFD methods that account for moving grids have been 
developed. Reference 1 1 presents a procedure for modeling configurations with moving 
rigid components using overset grids. In Ref. 12 a capability to model flows over flexible 
configurations with multi-block moving patched grids is presented. Based on these 
developments AIC computational accuracy can be increased by accounting for the 
effects moving grids. 

The AIC approach involves dividing the configuration into several panels based on 
stmctural details. Responses are then computed for several types of perturbations at 
different frequencies. For accurate computations each panel needs to be perturbed for 12 
modes, 9 strained and 3 rigid body modes [6], at about 5 frequencies. A typical space 
plane such as the X-37 [13] needs about 200 panels which leads to a total of about 12000 
cases for each Mach number. As a result, it requires enormous amount of computational 
resources. Such resources can be efficiently provided by single-image massively parallel 
supercomputers. 

NASA has facilitated development of large scale parallel computers since early 90’s [14]. 
Recently NASA in collaboration with Intel corporation and SGI has announced 
commissioning one of the world’s largest supercomputer that will be 10 times faster than 
the current ones[15]. As a result, computing CFD based AIC for aerospace design is 
becoming a feasible approach. 

This paper presents a procedure to efficiently compute AIC using the Navier-Stokes 
equations associated with deforming grids. 

Approach 

For a given configuration a suitable CFD grid is generated. From this the surface grid is 
extracted. Based on the CFD surface grid, AIC panels are defined. Each AIC panel 
is subjected to various modal motions to generate the aerodynamic coefficients 
and determine corresponding influence on other panels. Computation for each mode and 
frequency is assigned to a separate processor. 

Separate inputs are needed for each panel. An automated procedure that generates the 
panel input data is written in C++. Multiple cases are run using PBS (Portable Batch 
System) job control language[16] that has the ability to spawn cases to different 


processors. A flow chart of the process is given in Fig 1 . More details about the procedure 
will be given in the full paper. 
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Fig. 1 : Parallel AIC computing approach 


Results 


In order to test the process, a single panel rectangular wing with 5 modes is selected. 

The selected mode shapes are shown in Fig 2. For five modes AIC are computed at 5 
frequencies and 10 Mach numbers. AIC is generated using G03D, the 
Euler/Navier/Stokes streamwise upwind solver[17]. The grid deformations are computed 
at each CFD step using algebraic moving grids. Using the AIC data, a flutter boundary is 
computed. Figure 3 shows the good comparison of the computed flutter boundary with 
the wind tunnel test. 
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Fig 2 : Mode Shapes for a rectangular wing to compute AIC. 






Fig 3: Flutter Boundary of rectangular wing using AIC 

Next a multiple panel case is selected based on the simplified geometry of a typical 
orbital space plane[l 3]. CFD grid of size 600K grid points is generated for this 
configuration. The CFD surface grid that has 15000 nodes is shown in Fig 4. Using the 
C++ preprocessor AIC panels shown in Fig. 5 and corresponding input for G03D 
software are generated. 



Fig 4 : CFD Surface Grid for a Typical Space Plane 

The surface pressure coefficient distribution with a deformed panel in first bending mode 
is shown in Fig. 6. Each case using the G03D Navier -Stokes solver requires lOhrs of 
CPU time on a single SGI Origin 3000 processor. In the final paper AIC will be 
computed for all 200 panels by perturbing with 4 mode shapes at 5 frequencies. The 
flutter boundary will be computed. Scalability up to 4000 processors using the NASA 
Advanced Supercomputer Center[15], will be demonstrated. 


Fig 5: Panels of Space Plane for computing AIC. 



DEFORMED PANEL 

Fig 6 : Pressure Coefficient Distribution with Perturbed Panel 


Conclusions and Planned Effort 

A procedure to accurately generate AIC using the Navier-Stokes solver including grid 
deformation is presented. Preliminary results show good comparisons between 
experiment and computed flutter boundaries for a rectangular wing. A full wing body 
configuration of an orbital space plane is selected for demonstration on a large number of 
processors. In the final paper the AIC of full wing body configuration will be computed. 
The scalability of the procedure on supercomputer will be demonstrated. 
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